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Abstract. In this work we discuss on the basis of the 
instantaneous recycling approximation the possibility of 
determining a relation between matter distribution (lu¬ 
minous and dark) in elliptical galaxies and the observed 
metallicity gradients. 

A particular case is developed in detail yielding an 
analitycal expression of the above relation. 

Key words Galaxies: elliptical; Galaxies: chemical evo¬ 
lution; Galaxies: structure. 


1. Introduction 

Metallicity gradients in elliptical galaxies are inferred from 
the variation of metallic indices as a function of radius (see 
Danziger et al. 1993; Carollo et al. 1993 and references 
therein for previous work on the subject). The mechanism 
by which these metallicity variations could have arisen 
is still unclear. The models proposed are based on the 
processes developed during galaxy formation (Larson 1974 
and Carlberg 1984). 

Franx and Illingworth (1990) found that, for their sam¬ 
ple of 17 galaxies, the local (B-R) (U-R) colors are func¬ 
tions of the local escape velocity v esc for all of the galaxies. 
This result was confirmed later by Carollo et al. (1993) 
and Davies et al. (1993) using a more direct indicator 
of metallicity (M 52 ) than colors and a more appropriate 
calculation of v esc for each galaxy. This relation suggests 
that metallicity gradients could have arisen because of the 
different times of occurence of galactic winds in different 
galactic regions. In fact the galactic wind starts when the 
energy injected in the interstellar medium (ISM) by super¬ 
novae (of both type I and II) becomes equal to the binding 
energy of gas (Larson 1974). Moreover, it is reasonable to 
think that, when the galactic wind starts, the SFR breaks 
down. Hence, in the regions where v esc is low (i.e. where 
the local potential is shallow), the galactic wind will de¬ 
velops earlier and the gas is less processed than in the 


regions where v esc is higher, and the potential well in the 
external regions is lower than in the internal ones. 

This mechanism can produce metallicity gradients de¬ 
pending on the local binding energy of gas and therefore 
on the matter distribution. On the basis of this idea Mar¬ 
tinelli et al. (1997) reproduced metallicity gradients in 
good agreement with observational data, by using a Jaffe 
distribution for luminous matter (Jaffe 1983) and a dark 
matter distributed on a diffuse halo. In their work they 
do not take into account the processes occurred during 
galaxy formation such as collapse or merging. They pro¬ 
posed a model starting with the total mass already present 
at the beginning. This is a good assumption for elliptical 
galaxies since they must have formed in a reasonably short 
timescale (Matteucci 1994, Renzini 1995, Ziegler 1997). 

The aim of this work is to exploit previous idea in the 
frame-work of analytical models (based on the IRA ap¬ 
proximation) for deriving a relation between metallicity 
distribution and matter distribution in elliptical galaxies 
within the effective radius R e . A general formula is de¬ 
rived in section 2. Then in section 3 a particular case is 
considered and discussed. 


2. Relation between matter and metallicity distri¬ 
bution 

In the following, we consider elliptical galaxies to have a 
spherical symmetry. We assume that at the initial time 
(t=0) these systems consist only of primordial gas and 
the star formation starts at the same time in all galactic 
regions. 

The galaxy is partitioned in zones having the shape 
of spheric wreath, and thickness Sr with — << 1. Each 
region is assumed to evolve under the assumtions of istan- 
taneous recycling, good mixing and IMF (</>(m)) constant 
in space and time. 
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2.1. Basic equations 


2.2. Relation 


Under previous assuptions the equation for the gas evolu¬ 
tion in each region is (Edmundus 1990) 

dg = —ads (1) 

where g is the gas mass and ads is the amount of gas that 
does not return into the interstellar medium when ds of 
the interstellar material forms into stars. 

The equation for the evolution metals is: 

d(gz) = (jp — z)ads (2) 


Chemical evolution in a given galactic region proceeds 
until galactic wind starts. Exactly eq. (4) holds until 
t < tewi where taw is the time of galactic wind oc¬ 
curence. When t > tew we must replace eq. (4) by = 0. 
Therefore, we can establish the following relation between 
metallicity in stars and taw from eq. (8) 

< z >* _ e avtGW - 1 - avt G w , g , 

p gavtaw — 1 ' ' 


where p is the yield i.e. the ratio between the mass of new 
ejected metals and the mass which remains locked up in 
low mass unevolving stars and remnants. 

Under the assumption that z « 1 equation (2) has 


the following solution 


* = - p,09 V 

If we consider a linear star formation, i.e. 

( 3 ) 

ds . . 

_= m n 

( 4 ) 

we obtain from equations (1), (2) and (4) 


2 = pant 

( 5 ) 

g = g(Q)exp(-avt) 

(6) 


Equation (5) yields the metallicity in the gas. Hovewer, 
we are interested to know the mean stellar metallicity. 
The mass-averaged metallicity of a composite stellar pop¬ 
ulation is defined following Pagel and Patchett (1975) as: 

l r s M 

z ^ ds ( 7 ) 

s{t) Jo 

It should be noted that, although the mass averaged 
metallicity represents the real average metallicity of the 
composite stellar population, what should be compared 
with the metallicities deduced from observations is the 
luminosity-averaged metallicity. In fact, the metallicity 
indices, measured from integrated spectra, represent the 
metallicity of the stellar population which predominates 
in the visual light. The transformation from metallic¬ 
ity indices to real abundances should then be performed 
through either theoretical or empirical calibrations relat¬ 
ing metallicity indices to real abundances. Galactic chem¬ 
ical evolution models (Arimoto and Yoshii, 1987; Mat- 
teucci and Tornambe, 1987) have shown that the two 
averaged metallicities are different, since metal poor gi¬ 
ants tend to predominate in the visual luminosity so that, 
in general, the luminosity-averaged metallicity is smaller 
than the mass-averaged one. However, the difference is 
negligible for galaxies with mass M > 10 11 Mq. 

From equations (4), (5), (6) and (7) we have 

< 2 >* _ e avt -1-ant 
p ~ e aut -l 


The values of tew depend on r. In fact in the region at 
galactocentric distance r, taw is given by the following 
condition: 

EthsN i^GW i r ) — R B g as (fGW > r ) (10) 

where Eth SN (taw , r) is the thermal energy of the gas in 
that region and Esgas(tGW ,r) is its binding energy due 
to the gravitational attraction. 

We want now calculate explicitly the energies appear¬ 
ing in eq. (10). 

2.2.1. E thsN (t G w, r) 

The total thermal energy in the gas at the time tQw is: 

rtGW 

EthsN(tGW,r) = / eth SN (tGW - x)RsN(x)dx (11) 

Jo 

where Rsn(x) is the SN rate (either type I or II) and 
eth SN (tsN) is the fraction of the initial blast wave energy 
which is trasferred by the SN into the ISM as thermal 
energy. A detailed description of the SN rate is available 
in Matteucci and Greggio (1986). 

In the IRA approximation we can write 

Rsn(x) = 4:Trr 2 SrFi/g(x, r) (12) 

where 

p 8 plOO 

F = A (j>(M)dM + I <j)(M)dM (13) 

JMBm J 8 

and A and Msm are defined in Matteucci and Greggio 
(1986). IRA approximation related to the SNI is question¬ 
able because the SNI progenitors have mass not very large 
and consequently their lifetime is not negligible. Hovewer, 
for small values of tow the contribution of SNI is neg¬ 
ligible with respect to SNII (see fig. 1 of Matteucci and 
Greggio 1986). 

Substituting (12) in (11) and by changing the variable 
of integration we obtain 

rtGW 

EthsN (tow, r) = 47rr 2 5rFv / e t h SN (y)g(tGW ~ V, r)dy( 14) 



A. Martinelli: Metallicity gradients as detectors of matter distribution in elliptical galaxies 


3 


2 . 2 . 2 . EsgasitGW, r) 


where r is related to r by the following equation 


We define Esgasit-GW, r) as the energy necessary to carry 
the gas in a given galactic region at galactocentric distance 
r to r + l, where l is a model parameter which could be 
determined by means of a dynamical approach. Hovewer, 
in the following dicussion, its value is not relevant. 

Let M(r ) be the total mass within the radius r. We 
have therefore: 

r +l M(r') 

E Bg as{tGW,r) = 4nr 2 6rGg(r,tGw) — TT~ dr ' ( 15 ) 

Jr ^ 


Substituting eq. (15) and (14) in eq. (10), by using eq. ( 6 ) 
and by deriving with respect to r we obtain finally 


M(r + l) M(r) 

(r + l) 2 r 2 


Fv d 
~G~dr 



e th SN ( x)e avx dx 


(16) 


This last equation provides a relation between matter dis¬ 
tribution and tew- On the other hand we can invert eq. 
(9) to have tew as a function of < 2 >*. By substituting 
this last relation in (16), we obtain the relation between 
matter distribution and metallicity distribution. 


3 T 

1 - —-r -7 =- ( 22 ) 

p e T - 1 v ’ 

Equations (21) and (22) can be used to obtain informa¬ 
tions on the matter distribution in the following way. We 
consider a given value of r. By eq. (22) we can calculate 
the corresponding value of r (this is possible because the 
function in eq. (22) is bijective). By substituting this value 
in eq. ( 21 ) we obtain the difference ^+ 1)2 ~ ■ 

We observe that the function /(r) = ^ 
is always positive. Hence we can conclude from equations 
(21) and (22) that M(r) increases less than r 2 and there¬ 
fore the total mass density pTot(r) = PLum(r) + PDark(r ) 
decreases more rapidly than r^ 1 . This result is true for 
all the elliptical galaxies for which metallicity gradient is 
negative (i.e. 7 > 0). Previous result is in agreement with 
a Jaffe space distribution for the luminous matter (Jaffe 
1983) 


PLum( x ) 


1 

r 2 (l +r) 2 


(23) 


3. Discussion and conclusions 


The integral in equation (16) depends on r via tew and 
e th SN ( x )- Hovewer if we adopt the Larson (1974) prescrip¬ 
tion for e t h SN (x) (eth SN ( x ) = °- le o = 0.1 10 51 erg) the 
dependence is only on tew and we obtain from eq. (16) 


M(r + l) M[r) 0.1 Fe„ r , 

1+Tif — = ~^a~ (17) 

where r = avtcw and the prime denotes the derivation 
with respect to r 

By deriving equation (9) with respect to r we obtain 


e r 


(e T — l) 2 < 2 >'» 
t — 1 + e~ T p 


(18) 


The dependence of < z >* on r can be reasonably chosen 
as 


where r is normalized to 1 at the radius containing half 
the total emitted light in the space. This radius equals 
r e /0.763. Therefore, in the range 0.1i? e < r < R e , pLum(r) 
decreases as r ~ 2 (i.e. more rapidly than r^ 1 ) 

To obtain quantitative information from equations 
( 21 ) and ( 22 ) we should calculate the parameters enter¬ 
ing in those formulas. We must refer to the abundances of 
any a—element because they are produced on short time 
scales and hence IRA approximation is good. On the other 
hand the values of 7 strongly change from galaxy to galaxy 
and they are also depending on the calibration formula to 
compute cc—element abundance (for example Mg) from 
the observed indices. 

In a future work we intend to use the above relation 
for deriving constraints on the presence of dark matter 
in elliptical galaxies within the effective radius, since the 
mass M(r) appearing in our formula includes both dark 
and luminous matter. 


< 2 >*= (3r 7 (19) 

in accordance with many observational data (Carollo et 
al. 1993, Davies et al. 1993) in the range 0.1 < < 1. 

The quantities (3 and 7 are determined experimentally and 
their values strongly change from galaxy to galaxy. 
Therefore equation (9) becomes 

< 2 >* _ 7 < 2 >* _ 7 e T - 1 - r 

p r p r e T — 1 

Substituting eq. (20) in (18) we obtain from eq. (17) 

M(r + l) M(r) Q.lFeoje T — 1 e T — 1 — r 

~^G 
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